DFTT 52/93 
October 1, 1993 



THE GENERALIZED SIMPLIFIED 
PARTON SHOWER MODEL 



R. UGOCCIONI, A. GIOVANNINI, S. LUPIA 

(V) Dipartimento di Fisica Teorica, Universita di Torino and 

INFN, Sezione di Torino, via Giuria 1, 10125 Torino, Italy 

ABSTRACT 

> 

q We explore the consequences of considering clans real physical objects in 

| 7 j the framework of a generalized version of the Simplified Parton Shower 

model for a single jet. We predict that the average number of clans at 
fixed energy grows linearly in rapidity and slowly decreases with energy 
in a fixed rapidity interval. 

(N 

]~~] Talk presented by R. Ugoccioni at the 

CO XXIII International Symposium on Multiparticle Dynamics 

^ Aspen, Colorado, USA, 12-17 September 1993. 

To be published in the Proceedings. 



Work supported in part by M.U.R.S.T. (Italy) under Grant 1992. 



THE GENERALIZED SIMPLIFIED 
PARTON SHOWER MODEL 



R. UGOCCIONI*, A. GIOVANNINI and S. LUPIA 

Dipartimento di Fisica Teorica, Universita di Torino and 
INFN, Sezione di Torino, via Giuria 1, 10125 Torino, Italy 



ABSTRACT 



We explore the consequences of considering clans real physical objects in the framework of 
a generalized version of the Simplified Parton Shower model for a single jet. We predict that 
the average number of clans at fixed energy grows linearly in rapidity and slowly decreases 
with energy in a fixed rapidity interval. 



1. The Generalized SPS model 

Clan analysis has been a very useful tool for interpreting experimental data in high 
energy collisions; through its use, two regularities were discovered: first, the average 
number of clans is approximately constant with energy in a fixed rapidity interval, and 
second, the average number of clans grows linearly with the rapidity interval at fixed energy. 
These results have not been explained so far by any theoretical calculation, although they are 
well reproduced by Monte Carlo simulations. In an attempt to remedy this lack of analytical 
understanding of experimental facts, we propose to generalize an existing model 1 , based on 
essentials of QCD, by assuming clans as real objects. This assumption should be contrasted 
with the statistical clan concept used in the original model. 

In the original Simplified Parton Shower (SPS) model 1 , we consider an initial parton 
of maximum allowed virtuality W which splits at virtuality Q into two partons of virtuality 
Qo and Q\. We require Q > Qo + Qi and Qmin = 1 GeV. We define the probability for 
a parton of virtuality W to split at Q, p(Q\W), which is normalized by a Sudakov form 
factor. The joint probability density V(QoQi\Q) for a parton of virtuality Q to split into two 
partons of virtuality Qo and Q\ is defined by V{QoQ\ \Q) = p{Qo\Q)p{Q\\Q)K{Q), where 
K(Q) is a normalization factor. This general scheme is valid for any splitting function 
p(Q\W) which is factorizable in terms of its variables Q and W. However, for numerical 
simulations we choose: 
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A being the only free parameter of the model. This form of p(Q\W) was motivated in 
by our request of simplicity in the structure of the model. Notice that eq. (1) is exactly the 
virtuality dependence of the standard QCD kernel controlling gluon emission by a parton. 
For describing the rapidity structure of the model, we use the singular part of the QCD 
kernel controlling gluon branching: 

P(zi)dzi oc ( — + J dzi . (2) 

\Zi 1 - Zi ) 

Here z; is the energy fraction of the i th produced parton (i = 0, 1). 

We propose now to incorporate in the original SPS model which we have just 
sketched the concept of clan as follows: we decide to pay attention for each event to the 
ancestor which, splitting n times, gives rise to n subprocesses (one at each splitting) and 
we identify them with clans. Therefore in this model for a single event the concept of 
clan is no more as it has been in 2 , i.e., only a statistical one. In the present picture the 
number of clans in each event coincides with the number of splittings of the ancestor, i.e., 
with the number of steps in the cascade. It will be shown that this simple consideration is 
enough to predict that the distribution of clans in full phase space is Poissonian; we stress 
that this result can be obtained a priori in the generalized version of the model, differently 
from what has been done in 2 where the independent production of clans was an "ansatz" 
introduced a posteriori in order to explain the occurrence of NB regularity. 

The generalization of SPS to a genuine Markov process proceeds by assuming 
that the splitting function of the first step is the same at any step. This fact implies 
that our kinematical limits are extended from a triangle (Qo + Qi < Q) to a square 
(Qo < Q-, Q\ < Q), so that subsequent steps are not correlated. We call this new version 
of the SPS model Generalized Simplified Parton Shower (GSPS) model. It follows that the 
generating function of the number of clans in full phase space (fps) is given by: 

& n {z,W) = ze^ w ^\ (3) 

i.e., a shifted Poisson distribution. Accordingly, the average number of clans can easily be 
calculated by using eq. (1): 

iV(fps,WO = Ap(fps,WO+l = ^ P {Q\Q)dQ + \ = A\og(^^j+\. (4) 

In order to explain the clan scaling in rapidity, we extend our study to clan distri- 
butions in symmetric rapidity intervals Ay = [— y cut , y cut ]. We can relate the probability to 
have N' clans in the symmetric rapidity interval Ay, Pjy (z/cut, W), to the corresponding 
probability defined in full phase space, Pjv(fps, W), by the relation 

oo 

P N >(ycuuW)= Y, n(iV',y cut |iV,fps)P J v(fps,W), (5) 

JV=JV' 
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where Tl(N' , y cut \N, fps) is the conditional probability to have N' clans in Ay when one 
has N clans in full phase space. This new function contains all dynamical information 
on the rapidity structure of the production process. Since clans are not correlated in full 
phase space and the only allowed correlations are among particles within the same clan, 
the probability IT(iV', y cut |iV, fps) results to be a positive binomial distribution 

n(iV',y cut |iV,fps)= Q^^'a -*) N ~ N ', (6) 

where 7r(y cut , W) is the probability to produce a single clan in the rapidity interval Ay from 
the maximum virtuality W. Accordingly, it can be shown that the probability generating 
function for clans in symmetric rapidity intervals turns out to be 

9$L(z, W) = [7r(y cut , W)z + 1- 7r(y cut , W )]e x ^ w ^- l \ (7) 

with the average number given by Ap(y cut , W) = 7r(y cut , W)Ap(fps, W). 

Equation (7), being the sum of two Poissonian distributions (the first a shifted one), 
has a nice physical meaning: the two terms correspond to the probability of having the 
ancestor within or outside the given rapidity interval. In addition, it should be noticed that 
the void probability 3 in a given rapidity interval, Po (2/cut, W), is different from zero and is 
given by 

Po(y C ut, W) = [l- 7r(y cut , W )]e-^ W)X ^ W \ (8) 

This fact is far rich in consequences, which we propose to explore in a forthcoming paper. 
We also notice that, when y cut is very small and therefore 7r(y cut , W) tends to zero, the 
unshifted Poissonian dominates. Thus, it can be stated that in the smallest rapidity intervals 
the clan multiplicity is Poissonian to a good approximation. When y cut tends to full 
phase space and therefore 7r(y cut , W) tends to 1, we approach the exact, full phase space, 
shifted Poisson distribution. We find also that, when Ap(y cut , W) is sufficiently large, the 
shifted Poissonian dominates at large N (the tail of the distribution) while the unshifted one 
dominates at small N (the head of the distribution). This fact might have some consequences 
in interpreting the anomalies found in NB behavior for small N. Furthermore, the average 
number of clans in a symmetric rapidity interval, \y\ < y cut , can be calculated from eq. (7) 
and is given by the following formula: 

iV(y cu t, W) = 7r(y cut , W)N (fps, W) . (9) 
The next step is to calculate the explicit form of 7r(y cut , W), i.e., 

/2/cut 
p(y c ian, W)dy daD , (10) 
"2/cut 

where p(y c ian, W) is the probability to find a clan of rapidity y c i an generated at any step of 
the cascade initiated by an ancestor of maximum allowed virtuality W. The just mentioned 
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Figure 1. Average number of clans N(y cu t, W) 
as a function of rapidity interval y cut at different 
ancestor parton energies W in the GSPS model. 
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Figure 2. Average number of clans N(y cut , W) at 
partonic level as a function of rapidity interval y cut 
at different jet energies Ej for lowest energy jets 
from a sample of 100000 3 -jet events generated 
with JETSET 7.2 Parton Shower at ^/s = 1000 
GeV. Jet selection has been performed using the 
LUCLUS algorithm; partonic level has been ob- 
tained via GLPHD. 



probability can be expressed, within the GSPS, in terms of probability distributions calcu- 
lable from eq. (1) and (2). It results that the explicit form of 7r(y cut , W), valid in the central 
rapidity region, is given by 



7T(ycut,W0 = 



A 



Vcat 



A - 1 log W 

[(\og2W - y cut ) A - (\og2W + y cut ) A ] 



(ID 



+ 



2 A (A- \)(\ogW) A 



0< y cut < logW-log2 



By inserting eq. (11) and (4) into eq. (9), we determine the behavior of the average 
number of clans as a function of the rapidity interval and of the energy of the ancestor 
parton, which is shown for A = 2 in Fig. 1. A similar behavior occurs for other values of 
A > and our choice of A = 2 should be considered only indicative. As can be seen by 
inspection, the average number of clans grows linearly with rapidity at fixed energy and is 
a smoothly decreasing function of energy in a fixed rapidity interval. We point out that no 
one model to our knowledge has a so sharp analytical prediction on clan behavior in single 
jets. It is clear that our approximation consists in treating produced partons in a jet as pure 
gluons: our parton shower therefore should be considered a gluon shower. We notice that a 
similar behavior was seen by us 5 at hadronic level in single lowest energy jets disentangled 
by using the Luclus algorithm from a sample of 3-jet events generated with Jetset 7.2 
Parton Shower 6 . We remind that the lowest energy jets in 3-jet events are usually identified 
with gluon jets. In order to test our predictions in a region where data are not available, we 



4 



propose to reconstruct the parton level from the above mentioned sample of lowest energy 
jets in spite of the intrinsic limitations of Monte Carlo models and jet-finding algorithms. 
In view of NB regularity occurring in this sample 5 , we proceed by applying Generalized 
Local Parton Hadron Duality (GLPHD) prescription 1 . Results at the same jet energies 
seen in Fig. 1 are shown in Fig. 2: the predicted trend of the GSPS model is confirmed. 
The above predictions can be tested in the real world at LEP I and LEP II once a definite 
jet finding algorithm is universally accepted. Were our analytical predictions confirmed 
in experimental data, this fact would raise the question of the real existence of clans, i.e., 
groups of particles of common ancestor, in the production process and of their intimate 
physical nature. 
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